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Abstract 

We propose a QCD based many-body model for the nucleus where the strong cou- 
pling regime is controlled by a three body string force and the weak coupling regime is 
dominated by a pairing force. This model operates effectively with a quark-gluon La- 
grangian containing a pairing force from instantons and a baryonic string term which 
contains a confining potential. The unified model for weak and strong coupling regimes, 
is, however, only consistent at the border of perturbative QCD. The baryonic string 
force is necessary, as a stability and compressibility analysis shows, for the occurrence 
of the phases of nuclear matter. The model exhibits a quark deconfinement transi- 
tion and chiral restoration which are suggested by QCD and give qualitatively correct 
numerics. The effective model is shown to be isomorphic to the Nambu-Jona-Lasinio 
model and exhibits the correct chirality provided that the chiral fields are identified 
with the 2-particle strings, which are natural in a QCD framework. 

1 Introduction 

It has long been an important task for nuclear physicists to accommodate the most 
accepted fundamental quark interaction, QCD, in a field theoretical setting. This is 
because the nucleus constitutes the most important observational bound state system 
governed by the strong interaction. Furthermore, since the nucleons, the elements of 
the nuclear bound state, are believed to consist of quarks for which the rigorous strong 
interaction field theory applies, it is natural to try to describe the nucleus in terms of 
quarks (and gluons) beside nucleons. The present paper therefore aims at explaining 
nuclear properties and structure in terms of quarks and gluons. Still, one could be 
skeptical about using quarks with no direct observables in the description of nuclear 
physics that is very much based on experimental data. However, since the occurrence 
of strangeness excitations in hypernuclear physics needs to be explained by isospin or 
flavour quantum numbers that naturally come out of a quark model, there are good 
arguments for including quarks in nuclear physics. 
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There has been the obstacle of having a strong quark coupling constant above unity 
in the interesting range of the strong interaction which applies to the bound state and 
resonance phenomena. This strong coupling regime is characterized by low momentum 
transfer and corresponds to the region where quarks are believed to be confined within 
the hadrons, contrary to the high momentum transfer behavior where the coupling 
is small and quarks are believed to be free. A coupling constant above unity in the 
region of the strong interaction which is interesting for nuclear physics will basically 
render a perturbative description of the strong interaction field theory useless. We 
shall therefore use the trick of putting in quark clustering by hand, so that nucleons 
are treated as elementary particles, but quark structure is taken into account in the 
interactions by assuming that the mass of the nucleon depends on the medium through 
the local values of the fields a, tt. Once that is done we shall try to perturb the energy 
around the solution so obtained and else try to appeal to similar phenomena in other 
areas of physics (i.e. fluid dynamics) where solutions can be migrated to the nuclear 
system. 

In the beginning of the paper we shall introduce the quark description of the nuclear 
system and then propose a many-body Lagrangian for the nucleus. Then we shall go 
on with a description inspired on fluid dynamics and try to derive phases that could 
be relevant for nuclear physics. 

2 A quark model for nuclear structure 

In the following chapters we shall propose a general Lagrangian for the nuclear system 
in terms of quarks and gluons, and based on an interaction like that of gauge theories 
in the strong interaction version of quantum chrome dynamics, QCD. 

Basically the full Lagrangian is like a chiral model but we shall in the following list 
the various ingredients of the energy functional. 

2.1 The pairing force 

A starting point for a QCD model for the nucleus in terms of quarks and gluons could 
be a kind of a bag model in analogy to the well-known MIT model [1] for baryons and 
mesons but bearing in mind that the nucleons themselves are MIT bags. In the usual 
setting of the MIT model there is a strong influence of phenomenological aspects that 
we eventually will try to derive with field theoretical interaction terms. 

The generic form of the MIT bag model operates with a back pressure that ensures 
the quarks to be confined to the bag and a Coulomb force between the quarks. In 
the case, as here, where the entire nucleus is to be considered as an MIT bag with 
quarks moving within the boundaries of the bag, the bag pressure is again an infinite 
potential well that keeps the quarks confined to the nucleus. The quarks are, as usual, 
fermions with flavor and color charges. As for the interaction between the quarks, a 
crucial pairing force between the isospins ensure the total isospin to be zero. In the 
next subsection we introduce a baryonic string force that is responsible for nucleon 
formation and actually also for the bag pressure. 

The pairing force which has been derived from instanton contributions in QCD 
field theory [2] is able to account for important features of nuclear structure. For the 
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a) b) 

Figure 1: Baryonic string-configurations, (a) The angles of the triangle defined by a 3-quark 
cluster are all less than 2ir/3. (b)One of the angles of the triangle defined by a 3-quark 
cluster is greater than 2n/3. (c) A 6-quark cluster. This configuration is here included as a 
possibility but not used in the full picture. It will be discussed in a forthcoming paper. 




one-instanton contribution the pairing force is a truly 2-body force. Its expression is: 



Pi = gipi(x)(l + ij 5 )^i'(x)^j{x){l + i^ b )ip j :{x)e iljl e ij 

where anti-symmetrizes the isospin indices 

In the color zero sector, this effective interaction can be reduced to the Nambu- 
Jona-Lasinio form [3, 4], Pifoihh — Ctl7i) ■ (t^t!))- I n that case the missing terms, 
which constitute the difference between both expressions 

Pifo [(h + h!)(h + iil)(hh + ri • r 2 ) - (hi 2 - (f l7l 5 ) • (f 272 5 )) 

average out to zero in a mean field approach. 



2.2 The baryonic string force in the nuclear bag 



In this subsection we discuss the many-body string force term. We shall therefore 
introduce a three body string force to the Hamiltonian. This force replaces the MIT 
bag pressure and is supposed to be responsible for nucleon bag formation. The string 
force is proportional to the distance between the quarks being a very dominating force, 
so perturbation theory with the string force as a perturbation is not possible. On the 
other hand the strings will triplet-wise neutralize the color leaving us with only dipole 
or quadrupole moments. The string force is a qualified guess to a configuration - an 
ansatz. 

We shall now specify the baryonic string configuration. We have only limited infor- 
mation about nonperturbative string configurations. Lattice QCD gives, in the strong 
coupling limit, indication of color tube strings between two heavy quarks. To make a 
consistent picture of the pairing and baryonic forces we are staying within the border 
of perturbation theory, the strong coupling constant, o.qcd, being around 1. 

We can, however, resort to the dual string theories for help concerning string pa- 
rameters. Here, the string interaction is proportional to the minimal length L of the 
string: 

V st = -^-L = ^L (1) 
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where ol is the Regge intercept and a is the string tension. In the sequel, to avoid 
confusion with the sigma field, we write k = yfa. The total 3-body force is then, apart 
from color projection operators, 



The two possible [5] configurations between 3 quarks in SU(3) - QCD are shown in 
Fig. 1 when we require SU(3) gauge invariance. This means that quarks cluster in 
multiples of 3. Within our approximation we only consider two and three body forces 
arising in the color neutral configurations, Fig. 1 a) and b). 

Mesons are, as usual, strings between a quark-antiquark pair. 

2.3 The chiral aspects of the Lagrangian 

In this subsection we briefly mention relevant aspects of chiral symmetry, which we 
insist on having properly implemented. The pairing force is chirally symmetric. On 
the other hand, following the procedure of ref [8], the string interaction, which is not 
automatically chirally invariant, may be also reconciled with chiral symmetry. In our 
approach, the Dirac see is naturally responsible for the dynamical breaking of chiral 
symmetry, contrary to an artificial device such a Mexican hat. Our calculation shall 
be based on effective Lagrangians, either of the Nambu-Jona-Lasinio (NJL) type, or of 
sigma model type. We argue that both models are equivalent as far as the mean field 
description of bulk static properties of hadronic matter is concerned. However, they 
lead, in the RPA approximation, to manifestly distinct dynamics, in spite of being 
equally consistent with chiral symmetry, so that it really matters which model one 
should consider if, beyond bulk properties, we are interested on the surface properties 
or on the response of the system to an external probe. 



3 The effective action and the Dirac equations 



In this section, we describe the general structure of the effective action which is relevant 
for the present development and we try to motivate, on the basis of the Dirac equation 
of quark dynamics, the model which will be considered in the next section. The effective 
interaction is supposedly obtained by integrating over the gauge fields. 

Below we describe all the terms appearing in the effective Lagrangian, where the 
potential terms are coming out of a typical cluster expansion in the quark fields. First 
comes the kinetic term iipj^d^ip; next come all the 4-fermion terms making up the 
2-body force, 

which, following a common procedure, we, for simplicity, replace by the combination 
which is familiar from the NJL model [3], 



Finally, come the 6-quark terms contributing to the 3-body string force. In order to 
insure chiral invariance these three body forces can be understood as reductions of 
appropriate chiral invariant 4-body forces. 



V 3 = dx 1 dx 2 dx3ip(xi)ip(x 2 )'ip{x3)V s t'ip{x 3 )'ip(x2)ip(xi) 



(2) 



(3) 
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It is assumed that the string force is responsible for the transformation of the 
quark degrees of freedom into nucleon and a, fr meson degrees of freedom, in the 
strong coupling limit, at low energies. The a- nucleon coupling may be taken three 
times as large as the cr-quark coupling. However, as shown in Ref. [6, 7], the quark 
dynamics arising from the confining force has the consequence that the nucleon coupling 
constant decreases as the baryonic density increases. It will be argued that this effect 
is important for saturation of nuclear matter. 

In Ref. [6] the confinement of quarks into nucleons is described on the basis of the 
MIT bag model through the solution of an appropriate Dirac equation, 



(4) 



The a field is determined self-consistently, its source being the scalar density due to 
tp. If the MIT interaction Vmit [1] is replaced by a string force m(r) = nr, where k 2 
is the string tension, in a form which is consistent with chiral symmetry, the previous 
equation may be replaced by the following one [8] 



id t i> = -ia- VV> + (3(f n 1 m(r) + g) (a - ij 5 T ■ Tr)ip, 



(5) 



where a, n denote, respectively, the sigma meson and the pion field, and f n is the pion 
decay constant in the vacuum. Here, we have already factored out the three-body wave 
function into a product described by three independent equations. There are important 
differences between the models in Refs. [6] and [8] which explain the extra factors we 
introduce, namely (/~ 1 m(r) + g), and which are of interest for the present discussion, 
especially due to the fact that we insist on having chiral invariance. Moreover, the a 
field in our case couples to positive and negative energy states and therefore is nonzero 
in the vacuum. The relevant consequence of this, which will be explained later with 
respect to the phase analysis, is that we expect a much stronger saturation effect at 
high density of hadronic matter, and hence an enhancement of asymptotic freedom 
phenomena for quarks. 

We shall now check the consistency of our string configuration and the string po- 
tential m(r) by analyzing the spectrum of the Dirac operator: 



[a r p r + -a r (3k + /3<r(f n 1 m(r) + g)\ip = e(a)ip. 



where 



-i 

1 



1 
-1 



(6) 



(7) 



and k = ±1, ±2, • • ■. Here we have used spherical symmetry. We take, for r < lfm, 
g + f^ 1 m{r) = Kf~ 1 r, where k is the same as defined in connection with eq. (1). This 
gives us then, 
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F(r) 
G(r) 



(8) 



(e(a)-af- 1 K r)F + ^ + -G = 
dr r 

(s(a)+af- 1 K r)G-^- + -F = 
dr r 



(9) 
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Eliminating one of the variables we obtain a second order equation in the other which 
may be solved numerically. The equations will have the form 




= (10) 

and similarly for F. We have defined p = er. 

A few qualitative observations can now be made on the basis of this equation. This 
equation is singular for af~ l n = ±e 2 .This gives certain rules for the energy eigenvalue 
e(<r). It should be noticed that the string interaction is proportional to the a field. 

The QCD vacuum, which is studied in these calculations, is a very complex object 
not least because of the strong contributions from both positive and negative energy 
states. In our model it is described as a collection of nucleonic bags arising out of the 
solution of the Dirac equation, in all possible momentum states P such that |P| < 
A, where A is a regularizing cut off. In the framework of the independent particle 
approach which has been described, the mass of the nucleonic bag is approximately 
M(cr) = 3e(o"). The scalar density p$ and the energy of the vacuum are also important 
quantities and their calculation, which is left to another section, is very informative 
concerning this point. The vacuum value of the a field is set to /„-, i.e. a = f n if the 
Fermi momentum vanishes, pf = 0. 

The asymptotic solution for G and F fields in the Dirac equation can easily be 
found by approximating the equation to 

Kp dp \npdp I 
Thus the asymptotic solution for G and F will be 

G~exp(-g). 

Here we can easily read out the size of the nucleon bag to be which becomes 

approximately 0.8 fm, in vacuum. 



4 The model 

To be more specific, we present now in greater detail the basic scheme which supports 
our calculation. This framework is motivated by the previous discussion and is based 
on the assumption that quarks and gluons are the building blocks of nucleons. We 
assume that the interaction between quarks originates from an instanton force and 
from a confining string interaction, the latter being responsible for the clustering of 
quarks into nucleons and the former for the emergence of a sigma field which mediates 
the interaction between nucleons and is understood as a quark-antiquark disturbance 
of the vacuum. According to this picture, the well known sigma model Lagrangian 
density may be extended to read 

C = itp^d^tp - g s (a,7r)[(i)^)a + iitp^rip) • vf] + g v (tp^ip)uj^ 

-\ml(a 2 + vf 2 ) + \mluj^ + ^ad^a + ■ d,n) - ^u^, (11) 
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where ip denotes the nucleon field (not the quark field). This Lagrangian density de- 
scribes an assembly of nucleons, regarded as composite particles, which interact with 
a, tt and uj^ fields. The familiar Mexican hat does not appear in eq. (11). However, 
terms involving the w field are introduced, as in the Walecka model [9]. In a conven- 
tional model, with elementary particles, g s is a constant. In our model, which we call 
extended sigma model (ESM), the composite nature of the particles reflects itself in 
the fact that g s depends on the local values of the fields tt, a. We assume that 



1- Jl-^trV + T? 2 ) i . „ . . 

g.{*,*)=90\[ 2g 2 Q 2 ((j2 + 7f2) ~9o(l + ^° a ^ ( 12 ) 

where a, go are phenomenological parameters. For an extended system, tt = 0. The 
dynamically generated mass of the nucleon becomes then M(a) = g s (a)a. From eq. 
(12) it follows that 



M(a) = g Q a* 



l-Vl-4fl§aV „ . 1 2 



\ — — * 9oa{1 + 2 %a a ] - (13) 

We have circumvented the laborious problem of determining the mass M(a) of a 
nucleon subject to an external field a by making an ansatz. According to our as- 
sumption, M(cr) behaves as go a for small a, but increases faster for large a. We also 
have 

a 2 = —(1 - a 2 M 2 ). 
% 

The energy of an assembly of nucleons is 

£(a,u) = -r, ^ 2 + 9 2 * 2 -ri J2 W + \<<?V - 

PF<|P|<A \P\<PF 

or, after eliminating a in favor of M and "minimization" w.r.t. to, 

8(a) = -V E ^TW^ + ^M 2 {l-a 2 M 2 )V + 1 0^- A V. (14) 

p F <|P|<A Z9 1 6 m v n 

Here, V is the normalization volume, pp is the Fermi momentum, A is the regularizing 
cut-off momentum and r\ = 2Nf is the degeneracy. The natural assumption has been 
made that the wave function of the assembly of nucleons is a Slater determinant. 
Minimizing £(&), we obtain the generalized gap equation 

4?7 Y , 1 = 1 - 2a 2 M 2 . (15) 

The coupling constant depends on a, and, through a, on the baryonic density. The 
coupling constant increases essentially linearly with a 2 . The practical consequence 
of such a behavior is the emergence of a rapid increase of the energy with density, 
tantamount to the onset of a repulsive contribution. This behavior plays a decisive 
role in fixing the saturation density. 
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5 NJL versus sigma model 



In this section we discuss briefly the equivalence of the sigma model to the Nambu- 
Jona-Lasinio (NJL) model [3], as far as the mean field description of bulk properties 
of hadronic matter is concerned. The NJL model is defined by the Lagrangian density 



C = W-fd^il) + %[(^V) 2 + (V^srVO' 



(16) 



The sigma model is defined by eq. (11), provided we replace the coupling parameter 
g s (a,ir) by a constant. A regularizing momentum cut-off A is part of both models. 
Using standard procedures, the Lagrangian (16) leads to the Hamiltonian 



N 



G S 



N 



H NJL = ]T P k ■ a k + -± ]T 6(r k - r,)/3*/3i(l - ■ n). 



(17) 



k=l 



k,l=l 



The vacuum is described by a Slater determinant, |3>o)> created by the operators bp 
satisfying |P| < A, associated with plane wave negative energy eigenstates of the single 
particle Hamiltonian h = P ■ ot + (3M, and by the operators b P + , satisfying |P| < pf, 
so that is the Fermi momentum, associated with positive energy eigenfunctions of 
the same operator. The nucleon "constituent mass" M is a variational parameter. 
The energy expectation value £ = {^q\Hnjl\^q) reads 



V 



E — 

|P|=PF £P . 



T]V 

2^2 



A p4 
dP— 

PF £ P 



G S V 



V 
2vr 2 



A M p2 

dP 

Pf £ P 



(18) 



where ep = \/~P 2 + M 2 , V is the normalization volume, r) is the degeneracy and, 
for nuclear matter, Gs is 9 times bigger than the corresponding quark matter value. 
For quark matter, r] = 2N c Nf, while for nuclear matter, i] = 2Nf. The condition 
d£/dM=0 leads to the gap equation 



\P\=Pf £P 



V 



r]Gs_ 
2ir 2 



A p2 

dP— , 

sp 



(19) 



Pf 



which fixes M. For nuclear matter we fix A so that M = 939 MeV. For quark matter 
we fix A so that M = 313 MeV. If we take 

r 9 5o 2 

LrS = 5" ' 

this is essentially the same gap equation as (15), except that the quantity 1 — 2a 2 M 2 
in the rhs of (15) is here replaced by 1. 

Notice also that the minimum of £ given by eq.(18) is the same as the minimum of 
£{<j) given by eq.(14), if a is set equal to in the expression for M(a), eq.(13). 

An extended NJL (ENJL) model equivalent to (11) is easily obtained and reads 



+ ^-[{^f + (^ 75 rV) 2 ] 



£enjl 
K 

+y^[(V^) 2 + (#75^) 



Jl2 + %(^)(W)]- 



(20) 
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Table 1: Parameters of ESMI and ESMII 



9s 


9v 


A (MeV) 


a 


ESM I 5.81 


10.75 


477 


0.30 


ESM II 8.12 


13.16 


387 


0.25 



Table 2: Properties of nuclear matter according to ESMI and ESMII: binding energy, equilib 
density, incompressibility, /„-, chiral transition density and deconfinement density. 





E B 


Po 


M/M 


K 


fn 


Px 


Pde 




(MeV) 


(fm- 3 ) 




(MeV) 


(MeV) 


(fm" 3 ) 


(fm- 3 ) 


ESM I 


-15.75 


0.15 


0.7 


225 


150 


0.46 


0.687 


ESMII 


-15.75 


0.148 


0.58 


240 


108 


0.3 


0.36 



A similar model has been studied in Ref. [13]. 

6 The phases of nuclear matter and quark clus- 
tering 

In order to get a good understanding of the various phases nuclear matter can exhibit 
at different energies, it is useful to introduce some descriptors for the various nuclear 
processes. We usually describe quantum field interactions in terms of the momentum 
transfer of the strong interaction scattering involving quarks and gluons. One can as 
well picture the energy per nucleon as a function of pressure. 

There are two distinct transitions that are important in the nuclear phase diagram 
[10] when varying pressure or momentum transfer. One is the chiral restoration tran- 
sition, mostly determined by the instanton pairing force [2, 11]. Secondly there is the 
quark deconfinement transition when the nuclear "fluid" changes from "water" to be- 
coming "metallic fluid" . This is a consequence of the diminishing baryonic string force, 
predominantly due to the behaviour of the effective coupling constant. 

We wish now to calculate the energy per particle of nuclear matter as a function 
of the density. As pointed out before, our two main ingredients are the pairing force 
and the baryonic string force. We take these contributions separately since they apply 
mostly to different scales. The effect of the string force is to confine quarks within 
nucleons, while the pairing force comes from instantons. 

We have considered two extreme parametrizations of the ESM, denoted ESMI and 
ESMII. The corresponding results are shown in the tables and figures. The param- 
eters of the model, g s , g v and a, were chosen in order to reproduce the saturation 
properties of nuclear matter, namely, binding energy, saturation density and effective 
mass at saturation. The compressibility and are outputs of the calculation. We 
are now able to make the plots of the energy per nucleon as a function of density, 
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0.2 0.4 0.6 0.8 



P (fm- 3 ) 

Figure 2: Phase diagram of hadronic matter energy as a function of the density. The 
dotted line applies to a quark gas. The dash-dotted curve corresponds to the so-called 
NL3 parametrization of the Walecka model [15]. The full curve and the long dashed curve 
present, respectively, the results of ESMI and ESMII. The dashed curve corresponds to the 
quark- meson coupling model [6]. 

shown in Fig. 2. The dotted line applies to a quark gas described by the Lagrangian 
(16) with quarks. The coupling constant Gs = 2.14/A 2 was chosen so that /„- = 93 
MeV and A = 654 MeV. This cut-off reproduces the correct quark mass in vacuum, 
m q =313 MeV. We see that the quark equation of state exhibits no binding, although 
it shows a local minimum corresponding to chiral symmetry restoration. This occurs 
at around 0.2 fm~ 3 which corresponds to the generally accepted deconfinement tran- 
sition expected at a temperature of about 200 MeV in the vacuum. The full curve 
and the long-dashed curve represent the effect of quark clustering according to ESMI 
and ESMII, respectively. The cutoff A is such that chiral symmetry breaking of the 
vacuum reproduces the correct nucleon mass. The clear binding shown by these curves 
is due to the reduction in the kinetic energy arising from the clustering of 3 quarks into 
nucleons [12]. This latter behavior is crucial for stability. Our model, not only takes 
into account the effect of quark clustering, but also the repulsion between the nucleonic 
bags. To account for the interplay between these effects we have used eqs.(ll, 12, 14, 
15). In Fig. 2, chiral symmetry restoration is shown as a discontinuity of the second 
derivative of the ESMI and ESMII curves. Deconfinement occurs when these curves 
intercept the quark-EOS dotted curve. Although our calculation is indeed based on the 
ESM (eq.(ll)), we observe that it would equally be possible to use the ENJL model 
(eq.(20)), which also reflects the influence of the medium on the nucleon properties 
(see, for instance, ref. [13] or ref. [14], where a density dependent coupling is used, 
in contrast to our assumption of a coupling depending on the local value of the sigma 
field). However, we insist that the equivalence between the ESM (eqs. (11, 12)) and the 
ENJL model (eq.(20)) applies only to the description of bulk properties. These models 
differ in the description of surface properties and in the predicted dynamics, such as 
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Figure 3: Effective mass of the nucleoli and sigma field as function of the density for ESMI 
and ESMII. Chiral symmetry restoration occurs when these quantities vanish. 

the masses of scalar excitations, for which they lead to distinct dispersion relations. 

As has been argued in Ref. [6], the quark structure of the nucleon induces a 
mechanism for saturation by weakening, at high density, the attraction due to the 
sigma meson. We find that the contribution of this effect to the energy per particle 
is equivalent to a term essentially proportional to p 2 . In order to fit the high energy 
scattering nuclear experiment it is important that the slope in the high energy end 
(large pressure) is not too steep meaning a fluid of nuclear matter that is not too 
incompressible. 

In principle, one may expect to find two critical limits: the limit of chiral symmetry 
restoration, when the mass of the nucleon vanishes, and the deconfmement limit, when 
the baryonic string force becomes small and quarks get liberated. For ESMI, chiral 
symmetry restoration is depicted in Figs. 2, 3 and 4 at p/po ~ 3.1, while the transition 
from the confined phase to a deconfined phase is shown in Figs 2 and 4 at p/po « 4.5. 
For ESMII, chiral symmetry restoration is depicted at p/po ps 2.0, and the transition 
from the confined phase to a deconfined phase takes place at p/po ~ 2.4. We remark 
that the equation of state predicted by ESMI is consistent with the one predicted by 
the quark-meson coupling model [6], and is much softer than the one predicted by 
the so called NL3 parametrization of the Walecka model [15]. The equation of state 
predicted by ESMII is also rather soft. Our model is consistent with the occurrence of 
a superfluid exotic phase of hadronic matter between the critical points corresponding 
to chiral symmetry restoration and deconfmement. In particular, this phase might be 
composed by Bose particles, each particle being constituted by six quarks. Our model 
is also consistent with the occurrence of the so called CFL phase, above the second 
critical point. 

It may also be observed that the values of f n predicted by the ESMI and by the 
ESMII and shown in Table 2, respectively 150 and 108, are quite reasonable. 
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0.2 0.4 0.6 0.8 1 

P (fm" 3 ) 

Figure 4: Pressure versus density for ESMI and ESMII. Phase transitions are shown as 
discontinuities in the slopes of the curves. 

Finally, we wish to mention a recent attempt pursued by Chanfray and collaborators 
[16] to reconcile the Walecka model with chiral symmetry, although along different lines. 



7 Summary 

In the present article we have constructed a realistic quark model of nuclear structure 
phenomena which is inspired on the Quantum Chromo Dynamics (QCD) quantum held 
theory, being formulated with very few free parameters, namely A, a, go,g v . In its final 
form, our model resembles the conventional NJL model for nucleons with, however, 
consequential differences reflecting the influence of the medium on the properties of 
the nucleons. In order to meet important requirements of chiral symmetry arising in 
the description of nuclear phenomena at a sub-nuclear level, we have added to the 
starting Lagrangian which describes the dynamics at the quark scale, important terms 
coming from chiral fields that appear at a higher order of QCD, in a cluster expansion 
sense. Such a unified model of nuclear physics is consistent with important data about 
nuclear phases from low energy physics and astrophysics [17], such as chiral symmetry 
restoration and quark deconfinement. 
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Appendix: RPA treatment of the NJL and extened sigma 
models 

We argue that, although the extended NJL model and the extended sigma model 
proposed here are equivalent as far as the mean field description of bulk static properties 
of hadronic matter is concerned, they lead, in the RPA approximation, to manifestly 
distinct dynamics. 

For simplicity, we forget about the a dependence of the coupling coeficient g a . The 
Lagrangian (11) is chiral invariant. Therefore, one of the RPA modes comes at zero 
energy. This is the pion. Avoiding straightforward algebraic developments, we present 
the quadratic hamiltonian which describes the RPA modes having momentum, 



= ]T e k (cl T c KT + dl T d KT ) + m a (ala + b h° + bpb+ + b t 6 ) 



p F <|k|<A 



+ 



9 




k-r 



c k;r d_ k _ T )(a + aJ J ) 
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V*™ a v pF< | k | <A 
-c k , T d_ k ,_ T , (r T + ,(&+t + 6") + + 6+) + 2T T ° T ,(C + 6g)) 



Here, = \/k 2 + M 2 , , c, , c? are fermion and anti-fermion operators, a\a are 
scalar meson operators and 6°, b + , b^^ , b~ are pseudo-scalar meson operators. 
The quantities (Tj ,+T~ _,), i(Tj~ ,—T~ 2T^ r , are Pauli matrices, or, more precisely, 
they are tensor products of Pauli matrices in isospin times identity matrices in spin, 
the indices r being pairs of spin-isospin indices. Summation over repeated r, r' indices 
is understood. Terms responsible for the excitation of finite momentum RPA modes 
have been deliberately omitted. The quantities 

i^-k,-r> i^-k.-r^r" d_ k ,_ T ,c k)T 7? iT „ je{+,-0} 
behave as quasi-bosons. For instance 

[cL k _ r '»c k , r », 4 !T dl k _ T ,] T 3 T „ yil T^ T , ps <5 jfe ?? fc , 77+ = r/_ = 4, r? = 2, 

[d_ k _ T 'C kiT /,c kT dL k) _ T ] ps 4. 

The result requires replacing such quantities as dl k _ r ,cZ_ k) _ T '» and c kT c kjT » by 0, 
which is their ground state expectation value. Notice also that 

lcir>'CKT>>,ciJ_ K _ T ,]T T k T , = [dl % _ T J_^_ Tl ,J^J_ K _ T ]T^ T , = cl T dl K _ r ,T T k T „ 

As an example, we present now more explicitly the computation of charged pseu- 
doscalar modes. The following commutators are easily obtained 



4<? 



[n^,c KT d_^_ Tl T- rl ] 

[WW ,6ft] 
[W< 2 >,6bl 



:,t — k,— t r,r' 

-2e fc cg r a? 



.9 



4<7 



^' ' V / 2^F V "° 



-(b 



V (c£ T d f . -T+ - c k T d_ k , _ r T" 

CT p F <|fc|<A 



E (4,rrfLk,-r' 2 ?r' " Qk.-rd-k.-r'T; 



PF<|k|<A 



It follows that 



^ (2) , E ^(c^U,.^ + c k , T eZ_ k ,_ r ,T- r ,) - (6 - - 6 +t ) 

p F <|k|<A 

1 S ( 2Xk£fe + iSlrv) ( C lr d -*,-T> T ?,T> ~ CKrd-K-r'T-y) 

PF <|k|<A V V^rn a v / 



_%E k _ 

y/2m a v 



+ I (6 + bp) = 



VPF<|k|<A 

The result clearly follows from the gap equation, for 

9 



x k = - 



y/2m a V e k 
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This is the RPA pseudo-scalar mode and its eigenfrequency is 0. 

Similarly, we find that the eigenfrequency of the RPA scalar mode is close to m a . 
We present the computation of scalar modes. The following commutators are easily 
obtained 



[W (2) ,ir^k,-r] 
[W {2) ,C k , r d -k,-r] 

It follows that 



^/2m a V e k 



(a + a 



(«o + a o j ) > 



^2m a V Ek 



m a ap + 



9 



f^y E — (4,^ t -k,-r-Ck,rd-k,-r), 
^ lm ^ V PF <|k|<A £fc 



.9 



A/ 2^y E — (i^ t -k,-r-Ck,rrf-k,-r)- 
V^V* 7 PF <|k|<A £k 



p F <|k|<A 

~ E 2 ^(^fc4,r dt „fc 5 _ T - y fc C k,rd-k,-r) 
p F <|k|<A 

(V + Of l k l / t jt J A 



+t 



p F <|k|<A 



V- 4ff(A k + y k ) |k| _ +f _ +f 

+ 2^ / - = — (a +a ') +m (T (? ? a + ga ') 

p F <|k|<A e * 



= a 



E ( J kC kl r dt -k,-T + *kCk,rd-k -r) ~ ^0 + ■ 
\pF<|k|<A / 



The RPA equations read 



2£fcAk + ^ = — — ^ s A k 
2 £fc y k + ^±%M = _^y k 
^ 4g(A k + y k ) |k| 

p F <|k|<A 



^ 4ff(A k + y k ) |k| 



p F <|k|<A 



^2m a V £fc 



where X k , Y k ,£, 77 are the RPA amplitudes and f2 s the renormalized sigma mass. 
Let Gk = /0 2g — . The following dispersion relation is obtained 



2m a 



n 2 



m 



E si 



4£ fc 



CT p F <|k|<A 



n 2 - 4el 



k 



showing that f2 s rs m CT < 2M. Indeed, the l.h.s. of the dispersion equation becomes 
infinite for Q 2 = m 2 and for Q 2 = 4e 2 , changing sign at each place. The zero £l 2 
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lies slightly below m 2 . Between m 2 and 4M 2 , the l.h.s. of the dispersion equation is 
negative. Between 4M 2 + 4pp, and 4M 2 + A 2 we find a continuum of Landau damped 
modes. 

Next we consider the quadratic hamiltonian which, for the NJL model, describes 
the RPA modes having momentum, 

n( NJL = E £ k(. C l,T C KT + A,r d k,r) 
p F <k<A 



— ( C k,r d -k,-r ~ C Kr d -K-r)-r(cl,r' d -k',-r' ~ C k ',r'^-k' -r' 



p F <|k|,|k'|<A 

+ (4,T^-k',-T'^ryr' ~~ C k,r^-k,-T'^ T)T ') ( c k',r"^-k',-r"'^r",r"' ~~ C k' ,r" ^-k' ,-t"'T£, y. 



+ ( c k,r^-k,-r'^r,r' ~~ C k,r^-k,-r'^ r ,^ (4',T''^-k',-T'''^T",T''' - Ck' yd-^i-r'nT^, ^, 
+ 2 (4,r d -k ,-r' T r,r' ~ Ck,rrf-k,-r'^, T ') (4',r" d -k' -r"' T r",r"' ~ Ck', r "d-k' _ r '» T°„ r „, ) 

The following commutators are easily obtained 

— 16G S ( c k',r''^-k',-r'''^r'",r''' ~~ c k' ,r"d-k' ,-t'"T^, t „^ , 

p F <|k'|<A 

t^WJL' c k,rrf-k,-T'^ riT '] ~ -2efeCk ir (i-k,-r'^ 1 rir / 

— 16G S ^ ( £ Tc',r''^-k',-r'''^T",T''' _ Ck'jT"^-^",-^^ 1 ,-"^") • 

p F <|k'|<A 

Thus 

[WSk. E ^(ir^-k,-r^ T >+Ck,^-k,-r'T- T ,)] 
PF <|k|<A Z£fc 

16G 



E 



-7-^ E ( C k,r rf -k',-r' T rtr' " Pk^d-k.-r'^,) = 0. 

k J p F <|k|<A 



p F <|k'|<A £fc ' 

This is the pion mode. Notice that the NJL gap equation reads 

i- y. — = o 

p F <|k|<A £k 

and G s > 0. 

We present next the computation of scalar modes. The following commutators are 
easily obtained 



[ n{ NJLi C k,r d -k -rl ~ 2£ feC j c r d t _ k _ r - 16G S ^ ]T ^ ( C l ,r> d -k> ,-t> ~ C k',r'cU' -r') , 

p F <|fc'|<A 



£fc p F <|k'|<A £fc ' " ' " A "' 
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It follows that 



ni NJL> Z! ( X * C l,T d -k-T + y kC k ,r^-k,-r) 



p F <|k|<A 

-WG S (^k + lk)^- E ^(4',r^-k',-r'-Ck', T ^-k',-r') 

p F <|k|<A k p F <|k'|<A k ' 

= n s £ (*ki )T dL k _ T + YkC k , r d_ k _ r ) . 

p F <|k|<A 
The RPA equations read 

2e fc X jg - 16G S ^ 2 (X k , + y k ,)M = ft sXk , 

£fe PF <| k '|<A £k ' 



2e fc(^kC k>T (iL ki _ T - F k c k)T d_ k _ T ) 

p F <|k|<A 



2e fc y k - 16G S 



k'l 



p F <|k'|<A 



where -X" k ,Y k are the RPA amplitudes and S7 S the renormalized sigma mass. The 
following dispersion relation is obtained 



1 = -16G 8 J2 



k 2 4e fc 



PF <\k\<A £ k lL 

showing that Q s = 2M. Indeed, replacing S7 2 by 4M 2 in the l.h.s. of the dispersion 
equation we obtain an identity, in view of the NJL gap equation, 



-16G a ]T 



A: 2 



= -16G S E 



, . e? 4m 2 - 4el 

PF<|k|<A K ft 



k 2 Ae k 

PF<|k|<A K 



= 1. 



In conclusion, we have shown that both models generate pseudoscalar modes with zero 
mass. However, the sigma model generates a scalar mode with mass close to m a while 
the NJL model for nucleons generates a scalar mode with an unacceptable mass, equal 
to 2M. 



17 



